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Abstract 



We report on a non-perturbative determination of the parameters of the lattice Heavy Quark 
Effective Theory (HQET) Lagrangian and of the time component of the heavy-light axial- 

J> vector current with A^f = 2 flavors of massless dynamical quarks. The effective theory is 

k>< considered at the l/rrih order, and the heavy mass m-h covers a range from slightly above 

V^ the charm to beyond the beauty region. These HQET parameters are needed to compute, 

for example, the b-quark mass, the heavy-light spectrum and decay constants in the static 
approximation and to order l/mh in HQET. The determination of the parameters is done 
non-perturbatively. The computation reported in this paper uses the plaquette gauge action 
and two different static actions for the heavy quark described by HQET. For the light-quark 
action we choose non-perturbatively 0(a)-improved Wilson fermions. 

Key words: Lattice QCD; Heavy Quark Effective Theory 
PACS: 12.38.Gc; 12.39.Hg; 14.40.Nd 



1 Introduction 

Particle physics enters very exciting times with the start of cohecting data at the Large 
Hadron Colhder. Two ways are explored to probe New Physics (NP): either performing 
a direct search of new particles (e.g. at ATLAS and CMS) from the electroweak scale up 
to the TeV scale, or studying rare decays (e.g. at LHCb). The latter give rise to a very 
rich set of constraints on NP scenarios because they are either mediated by quantum loops, 
in which high energy particles circulate (that is the case, for instance, in flavor changing 
neutral currents), or by new current structures (in decays occuring at tree level). B-mesons 
offer a highly interesting and rich set of (rare) decay channels, such as 6 — )• 57, which arises 
in the standard model only from penguin diagrams, and hence puts strong bounds on NP 
scenarios. The analysis of inclusive decays such as -B — t- ^^7, in particular in the framework 
of a heavy quark expansion, strongly depends on the knowledge of the b-quark mass, m\y. In 
tests of the CKM mechanism, some tensions exist at the moment between sin 2/3, obtained 
from the golden mode B — t- J/^Ks, and the CKM matrix element V^b extracted from the 
B ^ Tu leptonic decay jljl. The theoretical input of the latter is the B-meson decay constant 
/b, whose uncertainty is ~ 10%. Currently, there is also a 3-0" discrepancy between V^^'^'^'^ 
and V^'^'^^'^ . Although it would be surprising if this difference were due to a significant 
underestimate of /s, reducing its error is an important task of lattice QCD. 

Lattice QCD will enable us to compute m\y and /b with a precision comparable to the 
one of forthcoming experimental measurements from high luminosity collisions. Nevertheless, 
a delicate issue for those extractions is how to get a satisfying control on the cut-off effects. 
Indeed, the Compton length of the b-quark is smaller than the typical finest lattice spacing 
of simulations in large volumes (used to compute hadronic quantities). Several strategies 
have been explored in the literature to circumvent this two-scale difficulty pHs], see [9] for a 
recent review. The ALPHA collaboration has proposed to use the framework of Heavy Quark 
Effective Theory (HQET) [lOpl with a non-perturbative determination of the couplings 12 



Its implementation at the first order in l/m,h has successfully been applied in the quenched 
approximation [13-16 . In this paper we will report on our effort to realize our HQET 



program for A'^f = 2 flavors of dynamical quarks, leaving the phenomenological results for 
forthcoming papers. Here, in particular, we present our determination of the couplings of 
the effective theory regularized on the lattice. 

We use the same notations as in [14] . The HQET Lagrangian density at the leading 
(static) order is given by 

^stat(x) = V^h(s)^oV'h(a;). (1.1) 

A bare quark mass rn^^^^ has to be added to the energy levels E^^''^^ computed with this 
Lagrangian to obtain the physical ones. For example, the mass of the B-meson in the static 
approximation is given by 

mB = i?^*^* + m^e • (1.2) 

At the classical level w-^are ^^ simply the (static approximation of the) b-quark mass, ?7ib, but 
in the quantized lattice formulation it has to further compensate a divergence, an inverse 



power of the lattice spacing. Including the l/?Tih terms, the HQET Lagrangian readaj 

^HQET(a;) = ^sta.tix) - UJi,i^O]iin{x) - UspinOspinix) , (1.3) 

Ounix) = Mx)'D^Mx), Ospin(x) = V^h(a^)cr-BV'h(x) . (1.4) 

At this order, two other unknown parameters appear in the Lagrangian, cjkjn and Wspin- Our 
normalization is such that the classical values of the coefficients are Wkm = ^spin = l/(2mh)- 
In order to compute the decay constant of a heavy-light meson, one needs the time component 
of the axial- vector heavy- light current Aq. At the lowest order of the effective theory, the 
current is form-identical to the relativistic one. At the l/jTih order, it is enough to add only 
one term to the static current (because we are only interested in zero-momentum correlation 
functions, see [TTJ ) 

^^Q^^(x) = Zl"^""^ [Ar\x) + c^l'^ 4\x)] . (1.5) 

In this work we present our non-perturbative determination of the parameters m-bare (the 
generalization of rn^^^^ to the l/rrih order), ln(Z^ ), Cj^ , 0;^!^, and Wspin at values of 
the lattice spacing relevant for the computation of hadronic observables. These parameters 
allow us to compute, for example, the spectrum of heavy-light mesons and heavy-light decay 
constants. As we explain in detail in the remainder of this paper, the basic idea is to match, 
in a small volume, a few observables expanded in the effective theory at finite lattice spacing 
to their non-perturbative continuum values determined in QCD at a given renormalization 
group invariant (RGI) heavy quark mass M. Collecting the observables in a vector $, the 
matching equation reads 

^HQET(^^ M, a) = $QC^(L, M, 0) , (1.6) 

where L is the space extent of the lattice and a — )• taken in QCD. Expanding the left hand 



side of eq. (1.6) at a given order of the inverse heavy quark mass defines a set of HQET 
parameters. We remind the reader that such an order- by-order treatment [12| is part of the 
very definition of HQET. 

The remainder of the text is organized as follows: in Sect. [2] we summarize the strategy 
of the computation, in Sect. [3] we give the details of its implementation, the final results 
can be found in Sect. [4] and Sect. [5] contains our conclusions. Some definitions are relegated 
to [Appendix A , in Appendix B| we explain how we tuned the parameters of the simulations 



and how we performed the renormalization of the QCD quantities, whereas Appendix C 
contains more technical information about the dynamical fermion runs. 

2 Strategy 

The computation reported here is done along the lines of [14]. For completeness, we repeat 
here the basic ingredients but refer the reader to this work for more detailed explanations. 
We start by the simulations of QCD in a small volume with space extent Li ~ 0.4 fm at four 
different values of the lattice spacing. We consider Nf = 2 dynamical light quarks that we 
tune to be massless and a quenched heavy (valence) quark that we simulate at nine different 



^The precise definitions of Do, D^ and ct-B can be found in 14 



mass values, such that the Ughtest mass is around the charm mass and the heaviest mass is 
above the b-quark mass. We compute five renormahzed observables that we extrapolate to 
the continuum 



^f"^^ (Li , M, 0) = lim $^Q^° {Li,M,a), i 

a— >0 



1, 



.5. 



(2.1) 



The definitions of these observables can be found in |14] and in Appendix A Here we just 



mention that ^i and $2 are finite volume versions of the heavy- light meson mass and the 
logarithm of the decay constant, respectively (up to kinematic factors). $3 is sensitive to 
the correction of the heavy- light current. Finally, ^4 and $5 are proportional to the kinetic 
and magnetic corrections, respectively. At the next-to-leading order of the effective theory 
(i.e. keeping the static and the 1/m.h terms) these observables can be expressed in terms of 
the five parameters discussed in Sect. [T| which we cast into a vector 



to 



rriha 



HZ, 



HQETx 



(1) 



(2.2) 



More precisely, we write the l/?n.h expansion of the observables in the following wa}r| 

$HQET^^^^^_ (2.3) 

The entries of the five-component vector i] and the five-by-five block-diagonal matrix (p are 
computed by performing a series of numerical simulations of HQET at fixed L = Li and for 



various lattice spacings a. A more explicit form of eq. (2.3) can be found in Appendix A 
As anticipated in the introduction, the matching condition that we impose is eq. (1.6) with 
L = Li: 

$^Q^'^(Li,M,a) = $Q^°(Li,M,0) . (2.4) 

Solving this equation defines the HQET parameters that we call w(M, a) 



(2.5) 
They are the bare couplings of the theory and, as such, can be determined by finite volume 



matching conditions. By imposing eq. (2.4), the parameters uji become functions of M, but 



this heavy quark mass dependence comes entirely from $Q*-^^. We then perform another 
set of simulations of the effective theory in a larger volume of space extent L2 = 2Li . The 



observables in this volume are then simply obtained by taking the continuum limit of eq. ( 2.3 ) 
in which we insert the parameters a)(M, a) computed in the previous step: 



$HQET(^^^^^Q) 



lim 

a-s>0 L 



ri{L2, a) + (p{L2,a) uj{M, a) 



(2.6) 



Our formulation of the theory and the non-perturbative determination of Co guarantees that 
all divergences, including those of order 1/a, l/a? are cancelled, and that the limit a — )• 
exists. Next, the parameters at larger lattice spacings (to be used in large volume) are 
obtained by inverting eq. (2.3) with L = L2, 



iu{M, a) = ip-\L2, a) ($«Qi^T(L2, M, 0) - r?(L2, 



a 



(2.7) 



^For example, in the case of the static-hght meson mass, r; is proportional to the static energy, and this 
equation is simply the finite volume version of me ~ E^^^^ + rribarc • 



Finally, in the last step we perform an interpolation (or, in one case, a slight extrapolation) 
in the inverse bare coupling /3 = G/^q and obtain uj{M,a) at exactly those values of the 
lattice spacing used in our large volume simulations |18|. 



3 Numerical application 

3.1 Continuum extrapolation of the QCD observables 

For the QCD simulation, we use the plaquette gauge action and non-perturbatively 0(o)- 



improved clover fermions 19 for Nf = 2 flavors of dynamical quarks with Schrodinger 
functional (SF) boundary conditions. We have four different values of the lattice spacing 
(< 0.02 fm), /3 varying in the range 6.16 — 6.64, and the number of lattice points per space 
direction being Li/a = 20,24,32,40. The physical volume is kept fixed by imposing for 
the Schrodinger functional coupling the value (^{Li/2) = 2.989, which corresponds to Li 



Ji 



0.4 fm. For each resolution Li/a we tune the dimensionless RGI heavy quark mass such that 

z = LiM G {4 , 6 , 7 , 9 , 11 , 13 , 15 , 18 , 21} . (3.1) 

With this choice, M varies approximately from 2 GeV to 10 GeV. More details about their 
renormalization and about the tuning of the bare coupling and quark masses can be found 
in [Appendix B For the run parameters we refer the reader to [Appendix C[ Finally, we also 



implement tree-level improvement, following exactly the procedure described in Appendix D 



of 14 



Our strategy for the continuum extrapolation differs somewhat from our previous work. 
The discretization effects can be important for our heaviest masses. In particular, for the 
simulations where Li/a < 24, we might have noticeable contributions of order (a/L)'^ with 
n > 2. We take advantage of the fact that we have various (and quite different) heavy quark 
masses: since the cut-off effects are smooth functions of a/Li and z, we perform a global fit 
of the form 

<^^^^{L,M,a) = ^^^^{L,M,0)\l + {a/Lif (A + Bz + Cz^)] , (3.2) 



using only the data points such that aM < 0.7, as motivated in 20,21 . Note that the two 



last terms in eq. (3.2) are proportional to (a/Li) x (aM) and (aM)'^. For each observable, 
<1>^ , the parameters of the fit are the nine different values $*°^ (L, M, 0) and Ai, Bi, Ci. 
As an illustration we show the results of the fit of <l>i and $2 in Fig. [T] We have checked 
that different fit ansatze (e.g. adding cubic lattice artefacts) give consistent results for both 
the central values and the errors, and that results are also compatible with the standard 
approach where the slope in (a/L)'^ is not constrained. 

3.2 Subtraction of the static part 

The effective theory is first simulated in the small volume of space extent Li, which is tuned to 
be the same as in QCE|^ We use five different lattice spacings, such that Li/a = 6, 8, 10, 12, 16 

^ In practice this tuning is done with a certain precision, which translates into a small error on the various 
observables. Since the static quantities are very precise this error is dominating for some of them and was 



taken into account as explained in Appendix B 
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Figure 1: On the top panels we show the continuum extrapolation of two QCD observables where we used 
a global fit to parametrize the mass dependence of the slopes in (a/L)^. The observable shown on the left is 
proportional to a finite-volume heavy-light meson mass and the one on the right to the logarithm of a finite- 
volume decay constant. On the lower panels we show the heavy quark mass dependence of these observables 
(z = LiM). 



(but for the continuum extrapolation we discard the coarsest point). The corresponding /3- 
values He in the range 5.26 — 5.96. For the static quark we use two different lattice actions, 
HYPl and HYP2, which are known to lead to small statistical errors [22]. Once again the 
light quarks are tuned to be massless, both in the valence and in the sea sectors. More 
details on the simulations can be found in the appendices. This set of simulations serves to 



compute the quantities 7/(Li,a) and ip{Li,a) in eq. (2.5). 

For i = 3,4, 5, the static contributions of $j are r/j. They have a well defined continuum 
limit and r/5 = 0. Thus we compute 



miLuO) 



limr]i{Li,a) 

a— >0 



3,4, 



(3.3) 



and show the continuum extrapolations in Fig. p) They are done linearly in (a/Li)"^, which 
is justified by two reasons: (i) The quantity 773 contains the time-component of the static 
axial current; we improve this current using the 1-loop value of ac^*** computed in [23]. This 
is sufficient since the improvement term has almost no numerical influence: even setting 
gives compatible results, (ii) The quantity r/4 is constructed from SF boundary-to- 



„stat 



boundary correlators, which are already 0(o)-improved in our setup, cf. Appendix C 
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Figure 2: Continuum extrapolation of rj3.{L\,a) and r]4,{Li,a), which are the static part of $3(Li, M, a) and 
'I>4(I/i, M, a) respectively. The two different colors represent the static action; blue for HYP2 and red for 
HYPl (slightly shifted to the right). 



e 











I 


0.2 








- 


0.15 




I 


i 


- 


0.1 


" i.-^^ 






- 


0.05 

n 











0.05 0.1 0.15 0.2 0.25 

l/z 




0.05 0.1 0.15 0.2 0.25 

l/z 



Figure 3: Mass behavior of $^ ^"^(Li, M, 0). As explained in the text, these quantities are expected to 
vanish in the static limit l/z = 0. 



Then it is interesting to define the l/w-h contribution to <1'3,4,5, 

$f/"^^(Li,M,0) = $f'^°(Li,M,0) -r/i(Li,0) , i = 3,4,5, 



and to study the mass behavior of these quantities, where $4 

-ll/ra) 



(1/m) 



(3.4) 
(Li,M,0) and 



$^^'"'''(Li, M, 0) are a pure kinetic and magnetic correction, respectively. The physical 



interpretation of $3 (Li,M, 0) is more subtle and involves the l/?n,h-correction to the 
time component of the axial current. Since they are expected to vanish like l/w-h at large 
mh) any strong deviation from a linear behavior could be interpreted as a contribution from 
higher-order corrections of the effective theor5Q One can see from Fig. [s] and Fig. |4] that 
our results for the heaviest masses are compatible with the expected (linear) leading order 
behavior in 1/ z. 

*Some deviations from linearity are expected since a 0(l/mh)-behavior always contains logarithmic modi- 
fications due to the renormalization of the effective theory. 



0.06 








I 


0.05 






i 


- 


Q 0.04 




i 




- 


o 


i 








aio0.03 


i 






- 


►Gh 


i 








0.02 


i 

i 






' 


0.01 











0.05 0.1 0.15 0.2 0.25 

1/z 

Figure 4: Same as Fig. [i] for tlie magnetic contribution <:E>*^'''"^(Li, Af,0) = $^°°(ii, Af, 0). 

3.3 Matching in Li 

With the set of simulations described in the previous section, we have computed ri{Li,a) 



and ip{Li,a). Thus the HQET parameters ijj(M,a) can be obtained from eq. (2.5), viz. 
However, in practice we split this equation in the following way 



Cbi{M,a) 



2 



.i(Li,a)(cI>QC°(Li,M,0) 



r]j{Li,a) 



+ 



5 

E 

i=3 



^~/{L,,a)[^f/"'\L,,M,0) 



(3.5) 



where we have used eq. (3.4). Whether one uses eq. (2.5) or eq. (3.5) to define iJj{M,a) only- 
affects the way we treat the lattice artefacts, but one expects a better precision in the second 
case. The reason is that the quantities r/3 and r/4 are by far the dominant part of <l>3 and 
$4 (and <I>5 vanishes in the static approximation). As explained in the previous subsection, 
they are static quantities which extrapolate to the continuum with O(a^) corrections. On 



the contrary, r^i and ??2 are divergent and have to be kept in the combination as in eq. (3.5). 
They are extrapolated linearly in a/Li. 

At this point we would like to comment about the separation of the different orders in 
the effective theory. The situation for the first two observables is different from the others, 
since 771 and r/2 do not have a continuum limit. However, the whole matching procedure 
can be carried out at static order as well. In that case, 0)3 = c]^ is just the improvement 
coefficient ac^^^^ which is approximated by perturbation theory. Moreover, the parameters 
(J4 and 0)5 are of order l/rrih, and therefore are set to zero. This setup defines the static 
approximation of the first two observables, ^f^^ and $|*^*, with cp = diag(Li, 1). Performing 



, 5) allows us to 



the matching only for these observables (instead of $^ with i = 1, 

determine the two parameters m^^^^ and ln{Z^^^) at static order. 

3.4 Evolution to a larger volume L = L2 

We then consider a set of simulations of HQET in which we use the same parameters as in 
the previous step (HQET in volume Li), but where we double the number of points in each 
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Figure 5: Continuum extrapolation of $i(L2,Af, a) and $2(^2, M, a) in the case where the 
done in the static approximation. We show the results for the nine different heavy quark masses 
discretizations HYPf,2. The conventions are the same as in Fig. IT] 



x10 

matching is 
and the two 



space-time direction. There we compute the quantities (p{L2, a) and r]{L2, a) with L2 = 2Li. 



We now have all the ingredients to compute $(L2,M, 0) using eq. (2.6), but for the reason 
given above, we use a slightly different version: 



#fQ^^(L2,M,0) = <' 



lim 

a-s>0 



Vi 



r]i{L2, a) + ^ VijiL2, a) ujj{M, a) 

5 
(L2, 0) + lim V] (pij{L2, a) ojj{M, a] 

a— >0 L ■^ — ' 



1,2 



i>3, 



(3.6) 



i=i 



Again, the continuum extrapolations of r/3 and r/4 are done linearly in (a/L2)^, while rys = 0. 



All other extrapolations are done linearly in a/L2 (see the discussions in Sect. 3.2 and 



Sect. 3.3). 



At the static order, the matching procedure at Li and the evolution to L2 can be done 
in the same way as with the five-component vector $. The result defines the quantities 
$^*^t(^2, M, 0) for i = 1,2, and their continuum extrapolation is shown in Fig. ^ When the 
matching is performed at the next-to- leading order, we have ^f^^{L2,M,0) = r]i{L2,0) for 
i = 3,4,5 and define the l/?TT,h-contributions as 

$f/'")(L2,M,0) = $f^''^(L2,M,0)-cDf'^t(^2,M,0). (3.7) 

Their continuum extrapolations are shown in Fig. [6] and Fig. [7j 

The parameters cj(M, a) are computed by the relation analogous to eq. (3.5) in the 
volume L2 



uJi{M,a) = J^v9r.i(L2,a)($fQ''^(L2,M,0)-r/,(L2, 
i=i 

+ ^v^r.i(L2,a)($f/-)(L2,M,0) 
i=3 



(3.8) 



In eq. (3.6) we have chosen to write the evolution of the observables to the volume L2 



in terms of the HQET couplings uji . Equivalently we can introduce a matrix of step scaling 



0.1 
0.08 
'CoO.06 



^H 



0.04 




0.02 0.04 0.06 0.08 

a/L2 



0.05 



S 0.04 



-Tf 



e 



0.03 
0.02 



• f — *— • • • 

, g — ,__, , • 

• • — »— • • • 

• m • — ■ • • 

, a a — m » • 

• ■ ■ — m » • 

• ■ • — ■ • • 



0.02 0.04 0.06 0.08 

a/L2 
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functions. In order to do that, one substitutes in eq. (3.6) the ljj by the matching eq. (2.5). 
One obtains an equation of the following form: 



$i(L2,M,0) = A^i(ii,M,0)+ lim Si(Li,a) 

5 

(1/m), 



+ hm VSi,(Li,a)$™(Li,M,0), i = l,2, 



J=3 



(3.9) 



$f/'"^(L2,M,0) = lim ySij(Li,M,a)$f/'^^(Li,M,0), i = 3,4,5. (3.10) 

' j=3 



The explicit definitions of D, S and S can be found in 14 . We have implemented the tree 



level improvement of these step-scaling functions in order to obtain a smoother approach 
to the continuum limit. Our results are obtained in this way, but tree- level improvement 
actually only has a small influence on our results after extrapolation. 



4 HQET parameters to be used in large volume simulations 

As we have already mentioned in the text and explained in detail in the appendices, the 
HQET parameters are obtained at five values of the bare coupling Qq, which are such that 
the renormalized coupling ^^(L2/4) is kept constant for the different ensembles. This is done 



by setting f3 to some precise values: /3 = 5.2638,5.4689,5.619,5.758,5.9631. In order to be 
able to use the HQET parameters in large volume simulations - and make phenomenological 
predictions - we interpolate (or, in one case, we extrapolate) them to ^ = 5.2, 5.3, 5.5, 5.7, 
using a quadratic polynomial in /3. We show this interpolation/extrapolation in Fig. ^ The 
results are reported in Table [T] for the action HYP2. All errors quoted in this table are 
statistical (obtained with a standard jackknife procedure), including errors coming from the 
QCD renormalization constants at finite lattice spacing. There is still an overall relative 
error contribution of about 0.9% from the quark mass renormalization in QCD, i.e., the 



factor h in eq. (B.5) of Appendix B, which relates the RGI mass to the SF running mass 
in the continuum limit. Since in practice this error will only become relevant when the 
z dependence of the HQET parameters is actually applied to interpolate to the B-meson 
scale and to extract physical quantities, it is enough to include it then. The complete set of 



results can be found on our website http://www-zeuthen.desy.de/alpha/, together with 



the relevant error-correlation matrices. From preliminary studies 18,24 , we know that the 
physical b-quark mass corresponds approximately to z = 13. For this value of z we display 
our results for ambare rescaled by Li/a for each value of the lattice spacing in Fig. [9J The 
parameter rribare absorbs the power divergences present in the binding energy of an heavy- 
light meson in HQET: a 1/a divergence at the static order, and a 1/a^ divergence at the 
l/?TT,h order. As one can see from the figure, m,bare is dominated by these power divergences. 
Clearly, in order to guarantee the existence of the continuum limit, these divergences have 
to be removed non-perturbatively, which is one of the benefits of our approach. 

Finally, we would like to comment on the size of 1/wi^ terms. In our setup, the fermion 
fields are periodic (in space) up to a phase. As a consequence, the static quantities depend on 
the choice of one angle, that we call Oq, and the quantities computed at the next-to- leading 
order of HQET depend on three angles: 6q,6i,92. So far in this work we have considered 



only what we call the standard choice of 6 angles 13,l4j, namely {Oq, 01,62) = (0.5,0.5, 1). 
When we change the values of Oi we change the set of observables, and thereby the matching 
conditions. The static quantities computed for different values of 6 are thus expected to differ 
by terms of order 1/mh- Once the l/rrih corrections are added, the difference should be of 
order l/rn^. In general we find no significant dependence for all the HQET parameters 
uji at the next-to-leading order of the effective theory, meaning that the l/m-^ terms are not 
visible within our statistical precision. As an illustration, we show the spread of our results 
for \a{Zf'^) and ln(Z^*^^'^) in Table |2| for the discretization HYP2, (3 = 5.3, and for z = 13. 
Due to statistical correlations, some of the errors on these 1/"^^ terms are significantly 
smaller than our errors for ln(Z||''^*) and ln(Z^ ) themselves, but still no significant 1/m,^ 
term is found. 

5 Conclusions 

We have reported on the first unquenched determination of a set of HQET parameters. 
While other approaches treat the dependence of the effective parameters on the renormalized 
coupling constant in a perturbative fashion, our approach is entirely non-perturbative in the 
QCD coupling constant and avoids the difficult-to-estimate uncertainties of perturbation 
theory for this case 17,25p Furthermore, our matching procedure takes into account non- 



^Ref. 125] computes the matching of the static-hght currents including 0(af ) and comments on a bad 
behavior of perturbation theory. In Section 3.3.2 of 17i, it is shown that different choices for the renor- 
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/3 


5.7 


5.5 


5.3 


5.2 




z = 11 


0.601(8) 


0.833(10) 


1.129(13) 


1.302(15) 


«"^bafc 


z = 13 


0.713(9) 


0.982(11) 


1.327(15) 


1.527(17) 




z = 15 


0.821(10) 


1.127(13) 


1.518(17) 


1.746(19) 




z = 11 


-0.209(4) 


-0.201(5) 


-0.197(5) 


-0.197(5) 


In(Zit^t) 


z = 13 


-0.191(4) 


-0.183(5) 


-0.179(5) 


-0.178(5) 




z = 15 


-0.176(4) 


-0.169(5) 


-0.165(5) 


-0.164(5) 




z = 11 


0.058(12) 


0.377(13) 


0.740(21) 


0.937(17) 


aJTlbaro 


z = 13 


0.236(12) 


0.582(14) 


0.985(17) 


1.207(18) 




z = 15 


0.395(13) 


0.769(15) 


1.212(18) 


1.459(20) 




z = 11 


-0.212(42) 


-0.194(37) 


-0.173(34) 


-0.162(33) 


ln(Z^Q^^) 


z = 13 


-0.181(40) 


-0.166(36) 


-0.148(32) 


-0.139(31) 




z = 15 


-0.155(41) 


-0.142(36) 


-0.127(32) 


-0.119(31) 




z = 11 


-1.00(14) 


-0.75(12) 


-0.63(10) 


-0.60(9) 


ci^Va 


z = 13 


-0.89(14) 


-0.68(12) 


-0.56(10) 


-0.54(9) 




z = 15 


-0.83(15) 


-0.63(12) 


-0.52(10) 


-0.50(9) 




z = 11 


0.778(13) 


0.608(11) 


0.485(8) 


0.441(7) 


Wkin/a 


z = 13 


0.682(13) 


0.533(10) 


0.425(8) 


0.386(7) 




z = 15 


0.609(13) 


0.476(11) 


0.380(8) 


0.345(7) 




z = 11 


1.626(59) 


1.284(49) 


1.041(39) 


0.956(35) 


Wspin/a 


z = 13 


1.404(51) 


1.109(42) 


0.899(34) 


0.825(30) 




z = 15 


1.236(47) 


0.976(39) 


0.791(31) 


0.727(28) 



Table 1: HQET parameters as a function of the bare coupling for the action HYP2 at z = 11, 13, 15. The 
central value 2 = 13 is close to the physical b-quark mass. The first two parameters, m^^^^ and \n{Z^'^*'), 
result from the matching at static order, while the remaining entries are the parameter set resulting from the 
matching of HQET at next-to-leading order. 



{Oi,e2 



Aln(ZJ^*^*) 



Aln(Z]JQET) 



(0,0.5) 



(0.5,1) 



(0,1) 



9o = 0.0 0.014(2) -0.046(47) -0.001(6) -0.012(14) 
% = 0.5 -0.048(49) -0.012(12) 

% = 1.0 -0.046(2) -0.055(60) 0.007(7) -0.008(18) 



ln(Zi*"*) 



Table 2: For each combination of 6-angles we compute the differences A\n{ZX 
ln(Zr*)|(e„,9i,e,)=(o.5,o.5,i) and Aln(Zj<^'='^) = ln(Z«'^'=^) - ln(Zj'^'=^)|(e„.e,,e,)^(o.5,o.5.i) with the HYP2- 
discretization at /3 = 5.3 and z — 13. As explained in the text, different values of 6 correspond to different 
matching conditions. The small ^-dependence observed at the static order is completely absorbed by the 
l/rrih corrections, meaning that the 1/m 
quantity. 



^ corrections are not visible within our statistical errors for that 
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Figure 8: Quadratic fits of HQET parameters in L2 as a function of /3. The blue points are the results of 
the numerical simulations and the magenta points are the interpolated (extrapolated) ones. In this plot we 
show the results for the HYP2 discretization and for the nine different values of the heavy quark mass. 
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Figure 9: The bare quark mass at the static order and its l/rrih correction, in units of Li, for z — 13 and 
the two discretizations HYP1,2. These quantities are expected to absorb the power divergences (1/a)", with 
n — 1 for the static order and n — 2 for the l/rrih correction. 



perturbatively the power-law divergences of the effective theory which can be numerically 
dangerous if simply addressed by relying on perturbation theory. We have shown that these 
divergences get absorbed into the effective parameters of the theory. The results presented 
here can be combined with hadronic matrix elements and energies computed in large volume 
simulations (some preliminary results have been reported in (18, 24, 261). We look forward to 



presenting our determination of the b-quark mass, of the B-meson spectrum and of heavy- 
light decay constants at the 1/w-h order of HQET in the two flavor theory, which will use 
the parameters computed in this work. 

Concerning the decay constant, the precision of the current matching is about 0.5% in 
the static approximation and 3% to first order in 1/mh. The latter puts a mild restriction on 
the available precision for the decay constant. It is reassuring to see that in the cases where 
we can check for contributions of l/m^ terms explicitly, these are considerably below our 
precision. This is in line with the strong hierarchy that we observe between different orders in 
the HQET expansion. Around the b-quark mass, the numerical values of <I>j(L2, M, 0) change 
from $1 « 20 = 0(L2mb) over $2 ~ 0.6 = 0(1) to <^[]/^^ « 0.006 - 0.06 = 0(l/mb), as 
one can see in Fig. [5]- Fig. [7} 

We finally note that our numerical computations have been carried out on apeNEXT 
computers, decommissioned by now. A future application of our matching strategy with 
three or more flavors can be expected to reach a further improved numerical precision on 
present or future hardware. 
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Appendix A Observables 

For completeness, we remind the reader of the definitions of the observables $ introduced 
in[l4], 

$QCD ^ j^^pP^ ln[-/A//A] , Ra, Ri, ^ln[/i/A;i]) , (A.l) 

where the quantities 

rP = -aoln[-/A(xo,^o)](,„=r/2,r=^), (A.2) 

Ra = In [/A(xo,6'i)//A(xo,6'2)](^^^y/2,T=L) ' (^•^) 



R, = iln 



/l (02)^1 (^2)3 



(T=L/2) 



(A.4) 



are built from the (renormalized) SF boundary-to-bulk correlation function /a of the pseu- 
doscalar channel and the boundary-to-boundary correlator /i and ki of the pseudoscalar and 
vector channel, respectively. The definitions are such that, at the 1/w-h order of the heavy 
quark expansion, the observables assume the following form 

cDf Q'^T = L(mbare + r^*^* + C^l^TsA + ^kinr>^^° + iO.^inT''^'^) , (A.5) 

^^^''^ = ln(Z^Q''^) + CA + 4'V,5A+a;kin*'^'" + ^spin^^P'^ (A.6) 

^HQET ^ i?^tat^^(l)^^^ + ^^.^^kin^^^^.^^^pin^ (^^^^ 

^HQET ^ Rf-^+uj^,^Rf-, (A.8) 

^^^''^ = C^spin^r, (A.9) 



which is just an explicit version of eq. (2.3). The precise definitions of the HQET quantities 



can be found in 14 . 



Appendix B Tuning of Li and renormalization in finite-volume QCD 

As explained in the main text, the basic element of our non-perturbative strategy to compute 
the HQET parameters consists in imposing matching conditions between a set of renormal- 
ized finite- volume observables in QCD, extrapolated to the continuum limit, and their coun- 
terparts in HQET, expanded up to order l/rrah- Since in this way the HQET parameters get 
determined by feeding results from non-perturbatively renormalized QCD into the effective 
theory, we summarize here how the renormalization in QCD (i.e., of the gauge coupling, the 
quark masses and the relevant composite fields) is performed. 

We work in a small volume of linear extent Li ~ 0.4 fm, where the SF setup [27|28| , with 
T = L and 9 = 0.5 as the periodicity angle of the sea quark fields, serves as our finite- volume 

14 



renormalization scheme. The physical volume stands in one-to-one correspondence to the SF 
gauge coupling (f'{L) running with the scale L 29,30 . We thus define Li by the condition 



f{Li/2) = 2.989 . (B.l) 

The known step scaling function of the SF coupling in two-flavour QCD t30j implies 

f{Li) = <7(2.989) = 4.484(48) , (B.2) 

and the associated exact value of Li in physical units could then be inferred from the results 



of 31 , but it is not relevant for the following. 

Our finite-volume QCD observables $^ , i = 1, ... ,5, are defined as suitable renor- 



malized combinations of SF correlation functions (see Appendix A) composed of a non- 
degenerate heavy-light valence quark doublet, where the light valence quark mass is chosen 
to be equal to the mass m\ of the mass degenerate dynamical sea quark doublet. The hop- 
ping parameter of the corresponding light valence (= sea) quark is denoted as k\ and that of 
the heavy valence quark by Kh- The $^ are universal and, in particular, their continuum 
limits exist, once they have been evaluated in numerical simulations along a line of constant 
physics specified by a series of bare parameters (L/a,/3, ki, Kh) such that the renormalized 
SF coupling and the light and heavy quark masses are kept fixed. 

According to eq. (B.l), /? is determined by requiring g^{Li/2) = 2.989 for given resolu- 



tions 2a/ Li. This peculiar value of the SF coupling was indicated by an initial simulation 
with Li/{2a) = 20, while for 10 < Li/(2a) < 16 additional simulations and interpolations 
in /3, based on the known dependence of the SF coupling and the sea quark mass on the 



bare parameters available from the data of 30 , were employed for fine-tuning to that tar^ 



get. Employing the non-perturbative /3-function of the SF coupling and estimates from our 
quenched calculation on the effect of propagating uncertainties in cf' (cf. Appendix D of flSJ), 
we can assess an uncertainty of about 0.05 in cf' to translate via Li into an uncertainty in the 
b-quark mass of at most 0.5%, which is negligible compared to the present direct uncertainty 
in the quark mass renormalization discussed below. 

In the quark sector, the sea and light valence quark masses are taken to be the same; in 
the numerical computation they are actually tuned to zero. The condition m\{Li/2) = is 
met by setting k\ to the critical hopping parameter, Kc, which is determined by the vanishing 



of the PCAC mass of the sea quark doublet, defined as in eq. (2.17) of 30 through the 



0(a)-improved axial current in the SF setup with boundary field "A" 29 . Again, Kc was 



estimated and partly fine-tuned on basis of the data published in [30], whereby for the 



improvement coefficient of the axial current, ca, 1-loop perturbation theory 32 and the 
non-perturbative estimates of [33] (after they had become available) were used. A slight 
mismatch of \-^m\{Li/2)\ < 0.05 of this condition is tolerable in practice. The resulting 



triples {Li/a, (3, Ki) are collected in the three leftmost columns of Table B.l Note that the 
mentioned 1-loop value for ca was only used in the preliminary determination of ki. All 
PCAC masses listed here are computed with the non-perturbative ca of [33|. 

It remains to fix the renormalized mass of the heavy valence quark to a sequence of 
values, fairly spanning a range from around the charm to beyond the bottom quark mass. 
To this end we choose the dimensionless variable z = LM, with M being the RGI mass of 
the heavy valence quark, because the latter is related via 

M = h{L) Zjnigo, L/a) (1 + &m(5'o) a"iq,h) mq,h + O(a^) (B.3) 
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L/a 



P 



Kl 



f{k) 



Zp{9o,t) 



Kh 



20 6.1569 0.1360536 2.989(36) 0.6065(9) -0.6633(12) 



24 6.2483 0.1359104 2.989(30) 0.5995(8) 



-0.6661(9) 



32 6.4574 0.1355210 2.989(35) 0.5941(10) -0.6674(23) 



40 6.6380 0.1351923 2.989(43) 0.5949(12) -0.6692(27) 



1.10443(17) 


4 


0.1327094 




6 


0.1309180 




7 


0.1299824 




9 


0.1280093 




11 


0.1258524 




13 


0.1234098 




15 


0.1204339 




18 


— 




21 


— 


1.10475(12) 


4 


0.1331966 




6 


0.1317649 




7 


0.1310257 




9 


0.1294907 




11 


0.1278628 




13 


0.1261106 




15 


0.1241815 




18 


0.1206988 




21 


0.1140810 


1.10455(17) 


4 


0.1335537 




6 


0.1325329 




7 


0.1320117 




9 


0.1309446 




11 


0.1298401 




13 


0.1286909 




15 


0.1274876 




18 


0.1255509 




21 


0.1233865 


1.10379(17) 


4 


0.1336432 




6 


0.1328462 




7 


0.1324413 




9 


0.1316178 




11 


0.1307738 




13 


0.1299065 




15 


0.1290126 




18 


0.1276125 




21 


0.1261232 



Table B.l: Bare parameters (L/a, /3, ki, Kh) used in the computation of the heavy-light QCD observables 
for L — Li. As explained in the text, they are fixed by the renormalization conditions (f{L/2) = 2.989, 
^ mi(I//2) ~ and LM = z for the SF coupling and the light PC AC and heavy RGI quark masses, respec- 
tively. The entering renormalization constants Zp and Z and the improvement coefficient fem were calculated 
in fail. 
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to the subtracted bare heavy quark mass, mq^hj and its hopping parameter, k^, in the 0(o) 
improved theory. Here, 

Z{go)ZA{go) „_ - 1 f 1 1 



^m(50, L/a) = ^§^f^ , amq,h = H ^ " ^ ) , L = L,/2 . (B.4) 



All ingredients of eqs. (B.3) and (B.4) are non-perturbatively known in two-flavor QCD: the 
axial current renormalization constant Za from [35|, and the renormalization factor Z and 
the improvement coefficient 6m from J34j. The scale dependent renormalization constant Zp 



for the specific /3- values in question was extracted in I34|, following exactly the definition 



of 36 . Zp, bra and Z are also listed in Table B.l In eq. (B.3), there also appears the factor 



^W = 4o = 1-521(14) , i, = l/L = 2/Li , (B.5) 

which represents the universal, regularization independent ratio of the RGI heavy quark 
mass to the running quark mass, m, in the SF scheme at the renormalization scale fi. 
h{Li/2) was evaluated by a reanalysis of available data on the non-perturbative quark mass 
renormalization in two-flavour QCD as published in l36i. 
Given the values 

LiM = z G {4, 6, 7, 9, 11, 13, 15, 18, 21} (B.6) 

of the dimensionless RGI heavy quark mass in Li and resolutions Li/a = 20,24,32,40, 



eqs. (B.3) - (B.5) can now straightforwardly be solved for the corresponding nine heavy 
valence quark hopping parameters Kh = Hii{z,go) that fix z to the numbers in eq. (B.6[^^ 



These hopping parameters and the associated z-values are collected in the two rightmost 
columns of Tab le IB. 11 

lists the bare parameters L/a, /3, k\, k^ (i = 1,...,9) of the numerical 



Table 



B.l 



simulations in L = Li, from which the heavy-light QCD observables <&i, i = 1, . . . , 5, are 
computed. These bare parameters are functions of the dimensionless variables (f'{L), z = 
LM and the resolution a/L, with well-defined continuum limits at fixed z, but we can also 
consider them as functions of the box size L, the RGI mass M of the heavy quark and the 
lattice spacing a. 



Let us still comment on the error budget arising from the above procedure of fixing z 34 
which has to be accounted for in any secondary quantity analyzed as a function of z. From 
the uncertainties on Za, Zp, Z, and fem quoted in the respective references i34l[35] (see also 



Table B.l ), one obtains by the standard rules of Gaussian error propagation an accumulated 



relative error on z in the range 0.38% < {di.z/ z) < 0.41% for all z-values and lattice resolu- 



tions in use here. The contribution from the universal continuum factor h{Li/2), eq. (B.5), 
represents with Ah/h = 0.92% the dominating source of uncertainty in the total error bud- 
get of Az/z = 1.01%. Note, however, that the error on this universal factor h has to be 
propagated into the QCD observables $, only after their extrapolations to the continuum 
limit. It is not included in the errors in Table IBTI 

A similar tuning procedure has to be performed for the HQET simulations necessary for 
the matching. Since the matching proceeds through renormalized quantities and therefore in 
the continuum limit, there is no need to use the same lattice resolutions on both the HQET 



* Owing to the sign and the order of magnitude of the non-perturbative values for bm in the /3-range relevant 
here, Kh(z,(7o) has no real solutions for arbitrarily high z- values. This implies that only for inverse lattice 
spacings Li/a = 24,32,40 hopping parameters Kh that achieve z — 18,21 can be found. 
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L/a /3 



Kl 



f{L) 



ami 



6 5.2638 0.135985 4.423(75) -0.01154(83) 

8 5.4689 0.136700 4.473(83) -0.00424(24) 

10 5.6190 0.136785 4.49(10) -0.00257(11) 

12 5.7580 0.136623 4.501(91) +0.00067(7) 

16 5.9631 0.136422 4.40(10) -0.00096(4) 

8* 5.4689 0.13564 4.873(99) +0.03189(18) 
12* 5.8120 0.136617 4.218(49) -0.00099(7) 

Table B.2: Bare parameters and results of the tuning to g^{L) = 4.484 for the HQET simulations entering 
the matching step. The additional lattices L/a — 8*, 12* are used to estimate and propagate a potential 
error, resulting from not meeting the line of constant physics condition exactly. 



and QCD sides. Larger lattice spacings can obviously be chosen in HQET compared to the 



relativistic case, as long as the choices are such that the condition in eq. (B.l) is fulfilled. 
The resulting bare parameters and values of g'^{L) and the bare PCAC sea quark mass, 

B we estimate 



defined as before, are collected in Table B.2 From the two runs at L/a 

^ ~ dz 



32=4 4g4 = 1.4(4), which we use to set a bound on the quark mass mi (which 
should be vanishing), such that its effect on the coupling is below the statistical error. The 
L/a = 12* label refers to an additional run used to propagate the error on g^ into the HQET 
observables. This is done by computing all HQET observables (at T = L and T = L/2) 
also at the bare parameters of 12* in order to estimate the variation of our primary HQET 
observables with respect to a variation of the renormalized coupling. This procedure neglects 
the lattice spacing dependence of the variation, which is justified for an error computation. 
In practice, the uncertainty on g^{L) is the dominating piece of the errors of 773, 774 shown 
in Fig. [2| while it can safely be neglected for all other observables within our present error 
budget. 



Appendix C Simulation details 

Our computations are in a natural way split into two parts, the generation of A'^f = 2 gauge 
field ensembles at the tuned parameters {L/a, (3,Ksca = ^i), and the subsequent computa- 
tion of all relevant correlation functions. The Sheikoleslami-Wohlert improvement coefficient 
Cswlffo) is set to its non-perturbative estimate |i9j for N{ = 2. In order to have an improved 
action in the SF we include boundary counterterms to cancel boundary-induced lattice arte- 
facts. The corresponding improvement coefficients, Ct((7o) ^ind Ct{go), are always set to their 
known 2-loop [37] and 1-loop 32 values, respectively. This guarantees that any boundary- 
to-boundary correlation function, such as /i, is 0(a)-improved. 

Ensemble generation: For simulating a doublet of mass degenerate, non-perturbatively 
improved dynamical Wilson fermions in the SF (^sea = 0.5), we use the algorithmic imple- 
mentation described in detail in |38| to produce the ensembles needed for the matching in the 
volume Li. Applying the step scaling technique in HQET to volume L2 = 2Li also requires 
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simulations at resolutions a/(2Li) while keeping all other parameters fixed. Furthermore, 
our strategy involves simulations at fixed time extent T = L as well as T = L/2. 

Since most of the simulations used for the HQET observables are fast and can usually be 
performed using several replica, we aimed for a total statistic of at least 8000 configurations 
in these ensembles. Only for our most expensive HQET ensembles with L2/a = 32 we did 
not reach this goal due to limited computing resources and thus restricted ourselves to have 
0(3000) configurations here. The QCD simulations have larger values of L/a (and thus 
smaller lattice spacings) compared to the HQET simulations. It is therefore more difficult to 
achieve high statistics for the QCD ensembles. Even more so since the gap in the spectrum 
of the Dirac operator, which allows to simulate at vanishing quark mass in the SF, decreases 
proportionally to the inverse time extent. Hence, for the production of ensembles used to 
measure QCD observables, our goal was just to reach a reasonable statistics, and thereby 
to obtain small and comparable errors in our final observables at the different resolutions. 
Thus, the ensemble size roughly increases from 0(500) at L/a = 20 to 0(1500) at L/a = 40. 

Using the notation introduced in |38j, we list the relevant algorithmic parameters and 



additional details in Table C.l For QCD the bare parameters {L/a, P,k\) are those of 



Table B.l, while for HQET we use those in Table B.2 which are not marked by a star. The 
molecular dynamics (MD) is characterized by specifying the trajectory length, the integrator, 
and the step size(s). For the trajectory length we choose r = 2 in MD units since we 
expect autocorrelation to be reduced |39] in that case. As integration scheme we always use 



multiple time scales with leap-frog integrator, also known as Sexton- Weingarten scheme 40 
The tunable algorithmic parameter po is introduced as a mass-preconditioning of the Dirac- 
operator a la Hasenbusch (41). The step sizes for the corresponding two pseudofermions 
are Stq and Sti, while for the gauge force we use 6To/6Tg = 4 throughout. (-/V^q) is the 
average number of conjugate-gradient iterations used to solve the symmetrically even-odd 
preconditioned Dirac equation during the trajectory, and -Pace is the acceptance rate of the 
simulation. Tmeas gives the MD time between configurations which have been stored on disk 
and used for measurements. In case we list more than one value of Tmeas, we have performed 
independent simulations with different measurement frequencies which have been chosen to 
be of the typical size of observed integrated autocorrelation times Tjnt. Explicitly, we show 
the average plaquette and PCAC mass of the production runs together with their estimated 



autocorrelation time in Table C.2[ We also list the results for the pseudoscalar boundary-to- 



boundary correlation function /i, which typically is the quantity with the largest integrated 
autocorrelation time among the different SF correlation functions. 

Measurements: since there is no explicit heavy quark mass contribution to correlation 
functions in HQET, we just need to specify ki = Kgea and the static quark action(s) in use 
to compute static-light correlation functions. The latter have been computed using the two 
static actions HYP1,2 described in |22| . For measurements in QCD we compute heavy-light 
observables in a partially quenched setup with k,\ = Kgea and heavy valence quark hop ping 

(i) I 

parameters Kvai,h = '^h ' ^ ~ 1, ... ,9. The latter slightly deviate from those in Table B.2 
because at the time we fixed the values of z the updated (7Q-dependence of Z^ fSS] entering 



eq. (B.3) was not yet at hand. This translates into a small mismatch (< 0.4%) of the z-values 
where we did the computation with respect to our target z-values. The QCD observables 
were interpolated to the target z- values to account for this. 
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sector 



Lja T/a [ 



(5to ' St\ J 



PO 



{N^D 



{N^ 



(e-^«) 



QCD[Li] 


20 


20 


40 


4 


0.0828 


51 


292 


87% 


0.9896(84) 


8 


i 


24 


24 


50 


5 


0.0755 


56 


346 


90% 


1.0006(45) 


10 




32 


32 


56 


4 


0.0651 


63 


442 


88% 


1.0009(48) 


10 




40 


40 


64 


5 


0.0450 


81 


533 


90% 


0.9986(42) 


4 




20 


10 


40 


4 


0.0828 


50 


145 


93% 


1.0012(29) 


8 




24 


12 


46 


5 


0.1140 


39 


166 


89% 


0.9991(56) 


10 




32 


16 


48 


4 


0.0977 


45 


211 


84% 


0.9990(83) 


8 




40 


20 


54 


4 


0.0870 


49 


257 


84% 


0.9983(80) 


6 


HQET[Li] 


6 


6 


30 


4 


0.150 


27 


106 


92% 


0.9989(8) 


10 


4 


8 


8 


30 


4 


0.130 


35 


137 


90% 


0.9983(13) 


10 




10 


10 


30 


5 


0.110 


32 


137 


89% 


1.0006(14) 


10 




12 


12 


32 


4 


0.100 


44 


191 


89% 


0.9990(14) 


10 




16 


16 


40 


4 


0.085 


50 


247 


91% 


1.0021(12) 


10 




6 


3 


24 


4 


0.243 


21 


47 


91% 


1.0004(12) 


8 




8 


4 


28 


4 


0.206 


24 


63 


90% 


1.0002(12) 


10 




10 


5 


30 


5 


0.182 


26 


77 


89% 


1.0016(15) 


8 




12 


6 


32 


4 


0.166 


28 


90 


89% 


1.0003(15) 


10 




16 


8 


40 


4 


0.141 


33 


118 


90% 


1.0003(13) 


10 


HQET[i2] 


12 


12 


52 


5 


0.100 


46 


256 


92% 


1.0018(32) 


6; 10 


i 


16 


16 


50 


4 


0.085 


52 


327 


90% 


0.9999(15) 


6; 10 




20 


20 


54 


4 


0.081 


54 


394 


88% 


1.0017(27) 


6 




24 


24 


50 


4 


0.070 


61 


441 


86% 


1.0003(23) 


6; 10 




32 


32 


64 


4 


0.063 


71 


747 


88% 


1.0059(78) 


8; 10 




12 


6 


40 


4 


0.166 


30 


114 


89% 


1.0010(14) 


6; 10 




16 


8 


40 


4 


0.141 


34 


144 


92% 


0.9996(25) 


6; 10 




20 


10 


50 


4 


0.100 


44 


177 


87% 


1.0006(13) 


6 




24 


12 


50 


4 


0.114 


40 


198 


85% 


0.9986(20) 


6; 10 




32 


16 


52 


4 


0.091 


48 


256 


86% 


0.9921(45) 


8; 10 



Table C.l: Algorithmic parameters of our production runs as explained in the text. 
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sector L/a T/a (plaquette) rint[plaq] 



(am) 



Tint [am] 



(/i) 



nntlfi 



QCD[Li] 


20 


20 


0.630066(13) 


3.6(6) 


+0.00055(13) 


8(2) 


0.4582(69) 


24(8) 


; 


24 


24 


0.636756(8) 


4.7(7) 


-0.000145(66) 


7(1) 


0.434(11) 


80(30) 




32 


32 


0.651034(3) 


4.1(5) 


+0.000146(32) 


5.4(7) 


0.4155(72) 


60(20) 




40 


40 


0.662409(2) 


3.9(5) 


+0.000034(17) 


2.7(3) 


0.3990(98) 


100(40) 




20 


10 


0.629675(14) 


4.8(6) 


+0.000182(57) 


4.4(5) 


0.9474(16) 


8(1) 




24 


12 


0.636422(11) 


5.0(7) 


-0.000400(56) 


5.4(7) 


0.9280(19) 


9(2) 




32 


16 


0.650779(7) 


5.0(9) 


+0.000001(36) 


4.1(5) 


0.8892(43) 


40(10) 




40 


20 


0.662203(4) 


5.6(9) 


-0.000056(22) 


3.8(5) 


0.8769(29) 


21(6) 


HQET[Li] 


6 


6 


0.546135(43) 


5.5(2) 


-0.00585(18) 


5.9(2) 


0.5041(11) 


5.9(2) 


; 


8 


8 


0.569268(25) 


5.1(3) 


-0.00339(13) 


5.9(3) 


0.5128(12) 


6.2(3) 




10 


10 


0.584467(14) 


4.6(2) 


-0.00260(9) 


5.7(3) 


0.5010(13) 


8.7(5) 




12 


12 


0.597377(10) 


5.1(2) 


+0.00040(6) 


5.4(2) 


0.4599(13) 


12.3(9) 




16 


16 


0.615017(5) 


4.9(2) 


-0.00107(4) 


5.7(3) 


0.4676(15) 


18(2) 




6 


3 


0.547894(68) 


4.6(2) 


-0.01686(11) 


4.7(2) 


1.1141(11) 


5.1(3) 




8 


4 


0.569319(34) 


4.6(2) 


-0.02096(10) 


5.4(2) 


1.0918(9) 


5.6(2) 




10 


5 


0.584064(23) 


4.7(3) 


-0.01167(7) 


4.5(2) 


1.0585(8) 


4.7(2) 




12 


6 


0.596924(14) 


5.0(3) 


-0.00338(5) 


5.6(3) 


1.0037(7) 


6.1(3) 




16 


8 


0.614569(7) 


4.5(2) 


-0.00226(3) 


5.7(3) 


0.9911(6) 


6.5(4) 


HQET[L2] 


12 


12 


0.546446(16) 


6.8(4) 


+0.00798(14) 


7.7(5) 


0.2527(17) 


19(2) 


i 


16 


16 


0.569646(8) 


5.2(3) 


+0.000457(65) 


4.8(3) 


0.3044(26) 


33(4) 




20 


20 


0.584826(6) 


4.9(4) 


-0.000900(55) 


4.3(3) 


0.3057(44) 


41(8) 




24 


24 


0.597708(3) 


4.7(3) 


+0.001342(32) 


6.2(4) 


0.2650(28) 


65(15) 




32 


32 


0.615288(5) 


3.6(6) 


-0.000841(61) 


4.6(8) 


0.280(18) 


130(60) 




12 


6 


0.547364(19) 


5.9(3) 


-0.003649(87) 


6.7(4) 


0.8076(9) 


7.4(4) 




16 


8 


0.569606(12) 


5.4(4) 


-0.002321(63) 


5.5(3) 


0.8047(10) 


6.6(5) 




20 


10 


0.584593(6) 


4.6(3) 


-0.001945(33) 


4.1(2) 


0.7840(10) 


10.0(9) 




24 


12 


0.597446(4) 


4.9(3) 


+0.000785(23) 


4.9(2) 


0.7120(9) 


14(1) 




32 


16 


0.615045(5) 


5.3(6) 


-0.000923(29) 


5.0(5) 


0.7322(25) 


29(6) 



Table C.2: Results for standard quantities measured during the production runs. 
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